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Abstract 



We consider evolutionary models for the population of short-period (j 10 hr) low-mass black-hole binaries (LMBHBs) and compare 
them with observations of soft X-ray transients (SXTs). We show that assuming strongly reduced magnetic braking (as suggested 
by us before for low-mass semidetached binaries) the calculated masses and effective temperatures of secondaries are encouragingly 
close to the observed masses and effective temperatures (as inferred from their spectra) of donor stars in short-period LMBHBs. 
Theoretical mass-transfer rates in SXTs are consistent with the observed ones if one assumes that accretion discs in these systems are 
truncated ("leaky"). We find that the population of short-period SXTs is formed mainly by systems which had unevolved or slightly 
evolved main-sequence donors (M? < 1 .2Mo) with a hydrogen abundance in the center > 0.35 at the Roche-lobe overflow (RLOF). 
Longer period (0.5 - 1) day)SXTs might descend from systems with initial donor masses of about 1 Mq and X,. < 0.35. Thus, one 
can explain the origin of short period LMBHB without invoking donors with cores almost totally depleted of hydrogen. Our models 
suggest that, unless the currently accepted empirical estimates of mass-loss rates by winds for massive 0-stars and Wolf-Rayet stars 
are significantly over-evaluated, a very high efficiency of common-envelope ejection is necessary to form short-period LMBHBs. 
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1. Introduction 

Soft X-ray transients (SXTs) are a sub-class of low-mass X- 
ray binaries most of which h arbour black holes (see, e.g., 
iRemiUard & McClintockl |2006|) . Their transient behaviour is 
commonly associated with the same thermal-viscous instability 
of accretion discs t hat drives ou tbursts of dwarf-novae (see, e.g., 
iDubus et al.ll200ll: lLasotall200U and references therein). In this 
article we will be interested in black-hole systems with orbital 
periods shorter than ~ 0.5 day. Some observational data on these 
systems are summarised in Table [T] All known low-mass black- 
hole binaries (LMBHBs) are transient. 

There are at least two major open questions concerning the 
origin and evolution of SXTs. First, the values of the parameters 
describing the common envelope phase, second the strength of 
the angular momentum loss through magnetic braking. 



1.1. The common-envelope phase 



has been challenged on the basis of computations which showed 
that envelopes of massive stars are very tightly bound to their 
cores. The "standard" equation for the variation of the orbital 
separation of components based on the balance between the 
binding ene rgy of the mass - losing star and the o rbital energy of 
the system (IWebbinklll984l: lde Kool et al. I ll987t) implies a ratio 
of final Of to initial a, separations of components equal to 
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As first suggested by iMcClintock & RemillardI (119861) and 
IdeKool etalTd 1987b the progenitors of LMBHBs may be rel- 
atively wide binaries (but still "close" in an evolutionary sense) 
composed of a massive primary (Mio > (25-40) Mq) and a low- 
mass companion (M20 $ 1 Mq). Such a binary avoids merging 
in the common envelope which is formed when a massive star 
overflows its critical lobe and survives the supernova explosion 
that produces the black hole. A black hole plus main-sequence 
star (henceforth, "bh-nms") binary is formed. As in cataclysmic 
variables, the further evolution of the system is controlled by 
the loss of angular momentum through gr avitational radiation 
and/or magnetically co upled stellar winds dde Kool et alJI 19871: 



, lly 

iPvlvser & SavoniieirT989.) . This evolutionary path for LMBHBs 



where acE is the common envelope ejection efficiency, A the 
parameter of the binding energy of the stellar envelope, M\ 
and Mi_c are initial mass of mass-losing star and the mass 
of its remnant, rn is the dimensionless radius of the star at 
the beginning of mass transfer, M2 is the mass of compan- 
ion. Formally, applying Eq. ([U to the estimate of the outcome 
of the common envelope stage one finds that a low-mass sec- 
ondary i s unable to unbind the envelope of a massive primary 
(Podsiad lowskieta l. 2003b; Kiel & Hurley 2006). Thus instead 
of formi ng a short- period binary, the comp onents will merge. A s 
found bv Podsiadlowski et al. (2003b.) and justham et al.l (|2006|) . 
producing a population of LMBHBs in a "standard" scenario 
requires the product oqe x /i to exceed ~ 0.1 which they con- 
sider as unrealistic. An alternative "anomalous magnetic braking 
scenario" dJustham et al . 2006) suggests that the progenitors of 
the donors of SXTs are intermediate-mass (> 2 M©) Ap/Bp-stars 
with anomalously high magnetic field strength. In this case, af- 
ter an initial high mass-loss rate stage of evolution, the secondary 
of the system turns into a low-mass star with a long evolution- 
ary lifetime. This scenario explains the short orbital periods of 
SXTs. However, in such a scenario the effective temperatures 
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Table 1. Estimates of spectral types, effective temperatures of donor-stars and mass ratios of components M2/MX in SXT. 



Object 


Po, 


Sp 




1 


Ref. Comments 




hour 











1 XTE J 1 1 1 8+480 (KV UMa) 4. 1 04 



2 GRO J0422+ 32 ( V5 1 8 Per) 5 .( 



3 GRS 1009-45 (MM Vel) 



6.840 



4 XTE J1650-500 7.680 

5 A0620-00 (V616 Mon) 7.752 



6 OS 2000+25 (QZ Vul) 



8.280 



7 XTEJ1859+226(V406 Vul) 9.120 

8 GRS 1124-68 (GUMus) 10.392 



9 H 1705-25 (V2107 Oph) 



12.504 



K7V-M0.5V 

K7V 

K5/K7V 

K5-M0 

K7-M0V 

mid to late K 

K5V-M1V 

M2 + 2V 
M2+^V 
MIV 
MIV 

later than G5V-K0V 
K7V-M0V 

K4 V 

later than a K3V, 
most likely 

between K5V and K7V 

K5 V 

K3V 



K5V 
K3-K6V 

consistent with K5V 



G5-K0 
K5V to K7V 
K3 -K5V 

K3-K4V 
K3/5V 

K7V 
K3V 
K5V 



0.083 



4700 ± 100 



4900 ± 100 



~ 0.008 
<0.1 

0.044 - 0.035 
116+°"™ 



0.313-0.076 



0.159-0.125 
0.1 



0.075 - 0.055 



0.053 - 0.035 



0.208-0.114 



< 0.053 



9 

10 

11 

12 

13 

14 

10 

15 

16 

10 

17 

18 



19 
20 

21 
10 
22 
23 

24 



10 
25 
26 
27 

28 
29 
10 

30 
31 
32 
10 



First IR-observartions 



B5 to K7 from {H - A:)o 
col. index 

possibly K6V 

next best 

match G5V and K2III 



evolved 

slightly evolved, 
but not a subgiant 
K4V is nearly identical, 
G5-K1 and K8-M0 also give 
very good correlations 

G5V fits best 

slightly evolved, 

K7 features not observed 



K3 to MO also give good 
correlations 



References: 1 -IWagner et aP ( 1200 ll) . 2 -|Zurita et alj l2QQ2|), 3 - iGelm o et aP JIOOl), 4 -i McClintQck et aJ.I i20Q3l), 5 - iTorres et all (l2004l) . 6 - 
[Mikolaiew.ska et al.72005), 7 - Gonzalez Hernandez et al. ( 2006), 8 - Casares et al. (2001), 9 - McClintock et al. (2001), 10 - Orosz (2003), 11 
-Casares et al. (1995), 12 - Harlaftis et al. ( 199^), 13 - Gelino & Harrison (2003). 14 - Revnolds et al. (2007) 15 - della Valle et al.iil997l) . 16- 
Filippenko et al. ( 1999), 17 - Orosz et al. ( 200 4|), 18 -Iproning et alj(l2007l). 19 - garrison et al. ( 2007) 20 - Shahbaz et al.1( ll999l) . 21 - 

^Gonzalez Hernandez et al. (2004), 22 - loannou et al.' ('2004^ 23 -'Harlaftis et aU lT99(^). 24 -.FiliDpenko et aL 41995,). 25 - 

Fihppenko & Chornock (2001), 26 - King et al. ( 1996b), 27 - Orosz et al. ( 199^). 28 - ICasares et al.l ( ll997h , 29 ^Shahbaz et ai] ( ll997h . 30 - 
Fili ppenkoetal.i(il997i) . 31 - iRemillard et all 11996) . 32 - iHarlaftis et al.i 01997 



of the descendants of the intermediate-mass stars significantly 
exceed the effecti ve temperatures of t he observed SXT donors. 
On the other hand lKiel & HurlevI (12006 ) found that a population 
of black holes accompanied by low-mass secondaries may be 
formed if the rate of winds from WR stars is reduced ad hoc. 

However, the estimates of both oce and A remain uncertain. 
The estimate of acE x /I is strongly influenced by the assump- 
tions about the mass-loss by stellar winds and by the uncertainty 
about the role of the internal thermodynamic energy in unbind- 
ing the envelope of the donor (see for the first suggestion of this 
source of energy in the context of formation of planetary nebulae 



i Lucvl (Il967h and e.g. iHan et al.1 (Il995l): iTauris & Dewil (l200lh: 
Sok er & Ha^ (l2003l) : iPodsiadlowski et all (l2003bl) : IWebbinkI 
(2007^ for further discussion of this issue applied to common 
envelopes). 



iDe Marco etall (l2003h carried out 3-D common envelope 
modelhng for an 1 .25 Mq AGB star engulfing 0. 1 and 0.2 Mq 
companions that took into account both rotation and the inter- 
action between the spiralling-in component and the donor star 
If the ratios af/oj obtained by De Marco et al. are inserted in 
Eq. ([TJ, they correspond to ckce x /I up to ^ 2. Regretfully, simi- 
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lar calculations are still absent for massive binaries and one has 
to rely on indirect methods for evaluation of uqe x 

The attempts to find a plau sible evolutionary sce nario 
for the pulsar PSRJ 2145-0750 dvan den Heuvell 11994 and 
the results of modelling of the population of binary pulsars 
(which also invokes high-mass stars ) both favou r gcE > < /I - 2 
(Portegi es Zwart & Yun^eisonlll998lrl iKalogeral (1 19991) in her 
study of evolutionary parameters of progenitors of donors in 
black-hole X-ray binaries found that the explanation of the ori- 
gin of the latter systems suggests qqe > 1, implying that sources 
other than orbital energy rele ase may be invoked in unbind- 
ing common envelopes. IHurlev et alj (1200 2) in their model for 
the total Galactic population of interacting binaries obtained 
a subpopulation of transient low-mass black-hole binaries for 
0.5 < Q-cE X ^ ^ 1-5. On the other hand using supernova 
rates and empirical es t imate s of the compact object merger rate 
lO'Shaughnessv et al.l (l2008l) constrained oqe x /I to be in the 
range 0.15-0.5. 

The issue of the sources of ene rgy that may increase q-ce was 
discussed by llben & Livid (11993! ). Referring the reader to the 
original paper, we mention only that, apart from recombination 
energy in ionization zones, Iben and Livio suggested, e.g., dy- 
namo generation of magnetic fields that may contribute to mat- 
ter ejection, enhanced nuclear burning due to injection of fresh 
matter into nuclear burning shells by circulation movements that 
develop in common envelopes and excitation of non-radial pul- 
sations that may drive mass loss. None of these mechanism has 
been explored as yet. 

All estimates of q-ce x A were obtained under different sets 
of assumptions on the evolution of massive stars that are con- 
sistent with our current knowledge of the stellar evolution. Thus 
it is still interesting to compare the predictions of various evo- 
lutionary models with observations. In this article we compare 
the Yungelson et al. (2006) model with the observed properties 
of secondary stars in SXTs (see also lYungelson & Lasotall2008h . 



1.2. Magnetic braking 



noted. 



Yet another probl em 
Ergma & Fedoroval (1199 
Ivanova & Kalogera (12006 
nism of AML hy_ low donor-mass 



e.g.. bv lKingetal.1 (|l996ah : 
iMenou et alf d 19991) and 
associated with the mecha- 
binaries. One finds that 



if, following IVerbunt & ZwaanI dl98ll) . one assumes that the 
braking law for single field stars Q oc f (S kumanichlll972h 
can be extrapolated over an order of magnitude in the rotational 
velocity v (from several 10 to several 100 km/s) to the case of 
close binary systems and if the spin-orbit coupling is efficient, 
the predicted mass-transfer rates for LMBHBs at orbital periods 
> 2hr are sufficiently high for these systems to have stable 
hot discs. However, such a population of stable and bright 
low-mas s black-hole X-ray binaries has not been observed. 
Also, the lSkumanichI dT972l) "law" is apparently in conflict with 
observational data on ro t ation velocities in youn g open clusters 
dCoUier CameronI l2002t lAndronov etal] l2003h . According to 
the latter study, the time-scale of rotational braking is two orders 
of m agnit ude longer than the one based on the Skumanich law. 
Also Ivan Paradiis (1996) noted that the values of M in SXTs 
are close to those expected if gravitational wave emission is the 
sole sink of angular momentum. Also for cataclysmic variables 
mass transfer rates predicted by the Skumanich law based AML 



disagree with observatio ns, see, e. g.. lHameurv et al.l dl988h and 
llvanova & Taaiiil d2003h . 

In our previous work dYungelson etalj|2006l henceforth. 
Paper I) we found that when the Verbunt & Zwaan AML mech- 
anism is allowed to operate after the systems become semi- 
detached, one obtains a large number of bright, steady LMBHBs 
that clearly are not observed in reality (also P. Charles, private 
communication). Therefore we suggested, in line with observa- 
tional evidence mentioned above, that in the semi-detached sys- 
tems with black-hole accretors, magnetic braking operates on a 
much reduced scale (as compared with the Verbunt & Zwaan 
prescription), or that it does not operate at all. As a test of this 
hypothesis, we computed a population of LMBHBs under the 
assumption that MSW is not operating once the RLOF occurs 
and have shown that in this case there remains in the Galaxy 
about the same number of such systems as in the case with ac- 
tive MSW (^ 10000) but all of them are transi ent, according t o 
the disc instability model (DIM) criterion of .Dubus et al.ld 19991) . 

In the present paper we extend considerations of the "no- 
MSW after RLOF" model proposed in Paper I to the case of 
ffCE x A < 2 and carry out a detailed comparison of the model 
with observations. 

2. The model 



For convenience we remind some basic information about our 
calculation of the LMBHB population. 

The model of the LMXB population is obtained in 
two steps: (i) modeling time-dependent formation of 
the population of bh+ms binaries, (ii) tracing the sub- 
sequent evolution of each system. The Galactic ensem- 
ble of bh-Hms binaries i s computed with the population 

synthesis code SEBA ( Porte gies Zwart & Verbunt '1996; 

Portegies Zwart & Yu ngelsonT 11998; Nelemans et al. 200l 
iNelemans et al.l 120041) using 250000 initial binaries with 
Mio > 25 Mq. The time- and position-dependent Galactic 
star formation hist ory in the code follows the model of 



' Though one cannot exclude that, because of the difference in the 
mass of progenitors of neutron stars and black holes, acE x A for them 
may be different. 



iBoissier & PrantzosI dl999l) ; for the inner 3 kpc of the Galaxy 
the star formation rate given in the latter stud y is doubled to 
mimic the Galactic bulge (see Figs. 1 and 2 of INelemans et al.l 
(2004)). The assumed binarit y rate is 50% (2/3 of the stars 
in binaries). The IMF follows iKroupa et alj (Il993l) . the initial 
distribution of semi-major axes of binaries (a) is flat in log a 
between contact and lO*" R©. A flat mass ratio distribution, and 
an initial distribution of eccentricities of orbits H(e) = 2e are 
assumed. 

For the common-envelope phase we used Eq. ([1]). We tested 
the combinations of common envelope ejection efficiency and 
stellar envelope binding energy parameters acE xA-2, 0.5, and 
0.1 (see below). 

Black hole progenitors have M[o = 25 - 100 M©. The rela- 
tion between MS masses of stars and pre-SN masses generated 
by SEBA agrees well with the one obtained by the SSE-code of 
Hurlev et al.. (.2000:) , despite differences in the treatment of stel- 
lar winds. The algorit hm for the formation of b lack holes follows 
the fall-back scenario (' Fryer & Kalogeral2001h with the assump- 
tion of a constant explosion energy of 10^" ergs, which is within 
the expected range. Nascent black ho les receive kicks at forma- 
tion that follow the iPaczvhskil d 19901) velocity distribution with 
cr,, - 300 km s scaled down with the ratio of the black hole 
mass to neutron star mass. 

In the next step of modeling, the population of bh+ms bi- 
naries born at different epochs is convolved with the grid of 
evolutionary tracks for low-mass components in the binaries 
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with different combinations of masses of components and post- 
circularization (initial) orbital periods (see Fig. |3]l. All tracks 
used in the paper were computed by an appropriate ly mo dified 
TWIN version (September 2003) of the Egg letonI (Il97ll) evo- 
lutionary code. As mentioned in ^1.21 the A ML via magnetic 
stella r wind was taken into account following lVerbunt & Zwaad 
([19811) : 



J = -0.5 X 10 



-2^-2^2 1^ 



M2RI 



(2) 



where M2 is mass of the secondary, R2 - its radius, ~ 0.1 
- its gyration radius, P - orbital period, / ~ 1 - a parameter 
derived from observations; it was set to 1 . Fo r momentum losses 
via gra vitational wave radiation the standard [Landau & Lifshitd 
( Il971h formula was applied. 

The evolution of each system was traced over a time span 
from formation to T - 13.5 Gyr or to the epoch when the mass 
ratio of the components of the system became q - M2/Mbh - 
0.02. At ^ i 0.02 the circularization radius of the accretion 
stream becomes larger than the outer radius of the accretion disc, 
resonance phenomena in the disc become important and it re- 
mains unexplored as yet how mass transfer then proceeds. The 
systems with q < 0.02 have f orb ~ 1 -5 hr, mass-transfer rates 
i 10"'" Mq yr"' and it remains to be observed if mass transfer 
occurs in them. 

For a more detailed description of the input parameters we 
refer the reader to Paper I. 



3. The population of progenitors 

As discussed above, indirect estimates of the possible range of 
the product of common env elope equation p arameters extend to 
acE X /l i 2. We note that iTauris & Dewil (1200 li) have shown 
that, depending on the definition of the core of the star and the 
treatment of the role of internal thermodynamic energy, the bind- 
ing energy parameter A may vary by two orders of magnitude 
(from 0.02 to 3.50 for the same 20 star at the tip of red giant 
branch) and therefore we consider q-ce x /I > 1 as an acceptable 
value. In Paper I we presented results of the modelling of the 
population of LMBHBs assuming a value of ckce x /I = 2. Here, 
we also discuss models with oqe x A-Q.\ and 0.5. 

A run of SEBA-code with cce x = 0.1 produced 
about 3400 zero-age bh+ms binaries formed in a Hubble time. 
However, all M2 i 1.5 Mq secondaries (i.e., stars subject to mag- 
netic braking) were paired with > 14 M© primaries, exceeding 
the largest estimate of bla ck hole mass in a known LMBHB 
(9.7 ± 0.6 for A0620-500, Froning et al.| [2007) - and often ex- 
ceeding the largest dynamically evaluate d mass of black h ole 
in binaries in the Galaxy (15.65+1.45Mq lOrosz etal.ll2007l) .R 
Since this model contradicts observations, we do not consider it 
further 

The model with ctqe x A-0.5 and AML via MSW implies 
the presence in the Galaxy of about 1700 semidetached bh-i-ms 
binaries with q > 0.02, 220 of which are bright and stable. 
Based on arguments presented in Paper I we conclude that, as in 
the acE X 't=2 case, in conflict with observations, up to several 
dozen persistent LMBHBs then would be observed in the Galaxy 
above the RXTE All Sky Monitor sensitivity limit, thus implying 
the need for a reduced MSW. The luminosities of these persis- 
tent LMBHBs are not high enough to correspond to the bright 
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Figure 1. Model distributions over separation of components in 
the population of precursors of bh+ms binaries with Mbh < 
12 M0, M20 < 1.2 M0, Porb < 1-2 day after circularisation of 
the orbits (upper panel) and in the ensemble of bh+ms binaries 
that produce LMXB (lower panel). Solid line - acE x ^^=2, dot- 
ted line - ffcE X /l=0.5. 



steady X-ray source s observed in elliptical galaxies (llrwi nl2006l: 
iSivakoffet al.ll2007h . 

Thus, we consider results for models in which AML via 
MSW does not act upon RLOF. Both for oce x A-2 and 0.5, 
almost all semidetached bh+ms binaries descend from systems 
with Mbh < 12 Mq and M2 < 1.2Mq. Table |2| compai-es the 
number of systems for two cases[3 The total numbers of sys- 
tems differ* by a factor ~ 2 and apparently neither contra- 
dicts observation- based estimates of several hundre d to several 
thousand objects ( Chen et alJll997l: [Romanilll998l) . though an 
ffCE X /l=0.5 model may be more attractive. 

In the Mio = (25 - 40) M© range to which most of the pre- 
cursors of black holes in bh+ms binaries belong, Mio » M20 
and the second term in the brackets in Eq. ([1]) is » 1. Then, 
roughly, after the common envelope stage one obtains the rela- 
tion: fl//flo °^ ffCE X A. Thus, the transformation law for compo- 
nent separation may be written, approximately, as f{a f)Aa j - 
C X /!(flo)Afl(), where hiao) and f{a /) are, respectively, the distri- 



^ The 34 hr bin ary X-1 in the starburst galax y IC 10 contains a > 
23 Mo black hole ( iSilverman & Filippenkdl2008h . 



^ Since we present one random realization of each model, all num- 
bers given are subject to Poisson noise. 
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Figure 2. Model distributions of initial masses of secondary 
components in the population of precursors of bh+ms binaries 
with Mbh < 12 Mq, M20 < 1.2 Mo, -Porb < 1-2 day after cir- 
cularisation of the orbits (upper panel) and in the ensemble of 
bh+ms binaries that produce LMXB (lower panel). Solid line - 
acE X ^-2, dotted line - uqe x /l=0.5. 

Table 2. Numbers of Galactic detached bh+ms systems with 
Mbh < 12 Mq, M20 < 1.2 Mq, Porb < 1-2 day formed in Hubble 
time, systems that reached RLOF, and current number of semide- 
tached LMBHB with q > 0.02. 





ttCE X 'i=2 


ttcE X A=0.5 


Total number of systems 


25685 


10910 


Systems that reached RLOF 


12150 


6150 


Systems currently alq > 0.02 


5080 


2980 



bution of systems over orbital separations prior to and after the 
common envelope stage. By virtue of this relation, functions / 
and h must be similar. This is illustrated in Figs. [T]and|2] 

For lower uqe, initially wider progenitor systems are sam- 
pled, as shown in the upper panel of Fig. [T] Initially too "close" 
systems merge, but they are replaced by initially wider systems. 
The rate of AML given by Eq. (|2]i is a function of orbital sep- 
aration (J oc fl""^) and as a result post-common-envelope sep- 
arations of components in progenitors of bh+ms systems that 
evolve into contact may descend only from a very narrow range 




Figure 3. Galactic zero-age population of low-mass black hole 
binaries (dots). The evolution of the systems that form LMBHB 
is determined by interpolation in the grid of evolutionary tracks 
with the border outlined by the polygon or analytically, if M20 < 
0.4 M0 and Pq < 0.8 day. Heavy dots mark initial parameters of 
the tracks listed in Table[3]and shown as examples in Figs.|4]-|6] 



of a, given a limited time-span from formation to the Hubble- or 
MS-evolution times. In this range of a the systems are distributed 
similarly irrespective of acE x A (Fig.[T] lower panel). 

Figure|2]shows that, while for acE x A=2 initial systems with 
low M2 dominate, "successful" progenitors of LMBHB have 
similar distributions over M2 both for acE x /l=0.5 and 2, since 
J oc M2/?2- Masses of black holes in the two oqe x A cases also 
have similar distributions, since for (25-40) M© stars that are 
typical of most of the progenitors of black holes in LMBHB, 
the masses of their He-cores do not change significantly in the 
hydrogen-shell burning stage that preceeds RLOF. As a com- 
bined effect of similar distributions of detached bh+ms progen- 
itors of LMBHB over Mbh M2 and a, scatter diagrams (e.g., 
M2-P) for populations of LMBHB are similar both foracE x A= 
0.5 and 2 and differ only in the number of systems per "unit area" 
(by factor ~ 2), since initial systems were sampled from different 
ranges in uq and our assumed initial distribution over separations 
is oc l/flo- For consistency with Paper I we further consider the 
model for acE x A=2. 

Figure |3] shows the distribution of zero-age ms+bh bi- 
naries in the initial-mass-of-the-donor M20 - initial (post- 
circularisation of the orbit) period Pq planeQif a bh+ms system 
reached contact, its further time-dependent behaviour was deter- 
mined by interpolaton in the grid of pre-computed evolutionary 
tracks. The borders of this grid are outlined in Fig. |3] The ini- 
tial systems with masses larger than the rhs side border of the 
outlined range or periods longer than the upper border of this 
region evolve to longer periods upon RLOF and have unstable 



Each dot in the plot represents several systems that have similar 
Mw, M20, Po but were bom at different epochs in the history of the 
Galaxy. Thus, some of the systems shown in the plot do not have time 
to evolve into contact. 
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Table 3. Parameters of tracks shown in the Figures (from left 
to right in each Figure. The columns list the initial mass of the 
star, initial orbital period of the system, period at RLOF, central 
hydrogen abundance at RLOF, the age of star at RLOF. In system 
1 the initial mass of the accretor is 12 M©, in systems 2-8 the 
initial mass of the accretor is 4Mq. 



No. 


Mo/Mo 


Pq, 




X, 


T, 






day 


day 




Gyr 


1 


1.0 


0.4 


0.316 


0.696 


0.07 


2 


1.0 


1.4 


0.375 


0.403 


4.33 


3 


1.0 


1.5 


0.395 


0.324 


5.30 


4 


1.0 


1.6 


0.420 


0.231 


6.33 


5 


1.1 


1.3 


0.425 


0.447 


2.59 


6 


1.1 


1.4 


0.446 


0.369 


3.22 


7 


1.1 


1.45 


0.460 


0.325 


3.51 


8 


1.0 


1.9 


1.772 


8 • 10-5 


9.17 



discs unless their mass is i 4 Mq. Evolution of the latter systems 
is illustrated by some evolutionary tracks in the Figures |4] and |6l 
Low-donor-mass systems, if their periods are not short enough, 
never evolve to contact. There is, however, a contribution to the 
population of LMBHBs from stars with M20 < 0.4 M© and ini- 
tial orbital period Pq < 0.8 day and these systems were evolved 
analytically. In binaries with M2() < 0.6 Mq and Pq > 0.8 day the 
donors do not fill their Roche lobes in the Hubble time. The ini- 
tial masses of the progenitors of the donors are typically < 1 
and this means that most of the donors have to be unevolved or 
slightly evolved at the instant of RLOF. 

4. Observational parameters of short-period 
LIVIBHBs 

4.1. Effective temperatures 

As noted by lJustham et al.l ( l2006l) . any formation scenario for 
LMBHBs has to explain the spectral types of black-hole low- 
mass companions. The determination of spectral types is a chal- 
lenging task since the emission of the cool star is contami- 
nated by the radiation from the accretion disc a nd the hot spot 
where the accretion stream hits the disc's edge dCharles & Cod 
I2OO6I) . Moreover, contamination by the disc may vary with 
time if the system did not reach quiescence. The published 
methods of spectral type determination vary in sophistication 
from naked eye estimates to using statistics after subtract- 
ing templa te spectra from Dopp ler-corrected averaged spectra 
(see, e.g., iHarlaftis et al.llT996l) : sometimes spectral types are 
inferred from colours, absolute magnitudes, SED and in many 
cases are subjective. For this reason, estimates of the donor's 
spectral type for the same system may differ by several subtypes 
(see Table [T] where we summarised the published spectral types 
for LMBHBs). 

The effective temperature of the donor is available in 
the literature for XTE Jl 118+480 (KV UMa) and A0620- 
00 (V616 Mon) only. They were derived as a by-product of 
abundance determinations that used synthetic spectra and x^- 
minimisation technique^ For other systems we were forced to 
apply the Sp -Teff relation for zero-age main-sequence stars. We 

5 iFroning et al.l |2007) cast dou bt upon the tempe rature determina- 
tion for A0620-00 by . Gonzalez Hernandez et aP | |2004.') . claiming that 
the latter authors used an insufficient set of spectral lines in their study 
and overestimated T^ff. 
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P, hr 

Figure 4. Model population vs. observational estimates of the 
ranges of effective temperatures of donors in SXTs (dots). 
Vertical lines mark the ranges of effective temperatures of donor- 
stars in observed SXTs corresponding to the ranges of the esti- 
mates of their spectral types (TablelTJ. Systems are annotated ac- 
cording to their number in Table [T] Large filled circles give Teft 
of donors derived from the fits to synthetic spectra. Heavy solid 
lines to the left and right show "limiting" tracks for a (l + 12)Mo, 
f = 0.4 day system in which MB does not operate after RLOF 
and a (l+4)Mo, Pq = 1.9 day system with MB operating after 
RLOF (see ^4. 11 for discussion). Thin solid lines are evolution- 
ary tracks for Mq = 1 .0 and 1 . 1 Mq donors with 4 M© accretors 
listed in Table|3] The Sp -T^ff relation used in the paper is shown 
at the right border of the coordinate box. 



used the relation given by Tokunaga in ICoxl (120001) . This relation 
is accepted by the authors of the catalogue of Hipparcos spec- 
troscopic standards. For A0-M6 spectra the approximate uncer- 
tainty of this scale (one standard deviation) is +100K. 

A caveat has to be entered concerning observers assigning 
"nonexistent" spectral subtypes to their objects. The modern MK 
classification system is devised in such a way that subsequent 
subtypes represent approximately equal differences in the spec- 
tra and some original decimal subdivisions were dropped. For 
instance, some subtypes between K5 and MO are absent: "K7 is 
considered a s half a subtype later than K5 and earlier than MO" 
(lKeenanlll985i) . For spectral types absent in the spectral classifi- 
cation (e.g., K3V) we applied a linear interpolation in log T^ff. 

In Fig.|4]we plot the distribution of our model population in 
the Porb - Tss plane and compare it with the T^ff of particular sys- 
tems. In the absence of objective criteria for the discrimination of 
reliable vs. non-reliable spectral type determinations, we plot for 
each system the complete range of the effective temperature for 
the range of spectral types assigned to it by different observers. 
Within the uncertainties of the spectral type determinations and 
conversion Sp -T^ff the model satisfactorily reproduces Teff of 
the donors in the LMBHBs with f orb ^ 9 hr. 

iPodsiadlowski et alJ (12003 a^ noticed that for Teff < 4500 K 
stellar models with grey atmosphere boundary conditions tend 
to overestimate the effective temperatures of stars, compared to 
models with more realistic non-grey atmospheres. Using un- 
evolved models of low-mass stars, Podsiadlowski et al. esti- 
mated that the correction in temperature may amount to about 
350 K. However, (i) they used a different stellar structure code, 
(ii) the correction may be different if the comparison is made 
for non-grey atmosphere models of out-of-thermal-equilibrium 
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mass-losing stars. However, models of the latter kind have not 
yet been computed and for the moment the difference between 
our models and the non-grey atmosphere models cannot be esti- 
mated. 

Some systems - GRS 1009-45, XTE 1650-500, A0620-00, 
and GS 2000-1-25 - are apparently located below the "populated" 
area. However, we restricted the initial periods of bh-nms bina- 
ries and masses of donors with the limits shown by the solid line 
in Fig. [3] so that interpolation between pre-computed tracks that 
upon RLOF evolve continuously to shorter periods and tracks 
that immediately go to longer periods (or change the direction 
of evolution in f oit) is avoided. As can be seen from the tracks 
plotted in Fig. |4]and from Table [3] (tracks 2 - 7), for a system 
with a given Mio and M20 the direction of evolution changes 
quite abruptly over a narrow range of initial APoib $0.1 day. 
For a given combination of Mio and M20 a "gap " between tracks 
evolving in different directions forms. But since there is a con- 
tinuity in the initial parameters of the systems, the "gap" in re- 
ality is filled. This is clearly shown by the tracks plotted in Fig. 
m Initial parameters of these additional tracks belong to a well 
populated area in the M20 - P[) diagram (Fig. [3). We did not 
pursue the goal of finding the precise parameters of the initial 
system(s) that may fit the parameters of a particular observed 
system and the precise borders of the progenitor space, since the 
accuracy to which the parameters of SXT are known and the un- 
certainty in the efficiency of magnetic braking do not justify this 
time-consuming and computationally expensive task. 

Nevertheless it is clear, at least qualitatively, that the origin 
of short-period LMBHBs may be explained within the paradigm 
of the strongly reduced magnetic braking in systems with donors 
overflowing Roche lobes. 

In Paper I we reduced the AML by magnetic braking (MB) 
to 0. This might be excessive and in reality some amount of MB 
can be still operating. (But we found that reducing the MB by 
a factor of 2 still leaves some ~ 100 bright steady sources.) We 
plot in Fig.mtwo "limiting " tracks: for (1-h12)M0, Porbo - 0.4 
day in which the donor is almost unevolved at RLOF and MB is 
absent after RLOF and for (1+4)M0, Poiho = 1.9 day in which 
the donor has Xc - 10 at RLOF and MB continues to operate 
(tracks 1 and 8 in Tab. [3]). Crudely, model populations with MB 
and without MB have to be located between these two limiting 
curves. Of course there will be a contribution from lower and 
higher mass systems, as we plotted only the IM© tracks for sim- 
plicity. Therefore adding some MB to our model will shift the 
population to the right, reaching a better agreement with obser- 
vations while still not producing stable luminous sources. From 
Fig. 6 in Paper I one can see that such an addition of MB mainly 
will influence the long-period systems. 

From Tableland Fig.|4]one can see that stars turn to longer 
periods if Xc < 0.35 at RLOF. For instance a I.IM© donor with 
^orbo = 1-4 (second to the right line) spends almost 10 Gyr in 
the RLOF-state, out of which during about 5 Gyr it evolves to 
longer periods. This provides the possibility of explaining SXTs 
with periods > 9-10 hr However, it is then necessary to compute 
a new grid of tracks for stars that evolve to longer periods and 
a very dense grid of tracks to cover the space between tracks 
evolving to shorter and longer Porb- This will be the topic of a 
dedicated paper. 




P, hr 

Figures. Mass transfer rates in model LMBHB as a function 
of their orbital periods (dots). Arrows mark upper limits to the 
estimates of mass-transfer rates in observed SXTs as given by 
Eq. Q. Crosses are estimates of mass-transfer rates in SXTs 
based on recurrence times. 



time. However, this approach has several weaknesses. First, re- 
currence times are known only for a few systems. Second, it is 
not sure that the rate calculated in this way represents the secu- 
lar value (this is the general drawback of mass-transfer rate es- 
timates) and third, it assumes that during the "refill", accretion 
onto black hole does not occur This last assumption is put in 
doubt both by observations (see e .g. Done et al. 20 07. and ref- 
erences therein) and models (see lLaso"tall2008l and references 
therein) which suggest that quiescent discs are truncated and 
therefore leaky. In such a case one can estimate the upper limit 
to the mass transfer rate which cannot be larger than the critical- 
for-stability accretion rate at the truncation radius (see Paper I 
for details). The actual mass transfer rate should be somewhere 
between the values estimated by the two methods. 

Figure|5]compares the model mass transfer rates with obser- 
vational estimates of M. We present two estimates of the latter. 
We show M estimates from recurrence times and mass accreted 
during the outburst and estimates of the upper limit of the accre- 
tion rate at the truncated disc inner edge. In the latter case, we 
get: 



2.5 ■ 10-^ [(1 + q^l^ (0.5 - 0.227 log^)] ' P^^fi-^^ Mq yr^S) 

In Eq. (O Pel is the orbital period in days, and /, < 0.48 is the 
fractional disc truncation radius. We used the revised version of 
the critical accretion rate (iLasota et al.ll2008l) : 



M'rit = 2.64 X 10'^ Rf^^ Mi-" **5 g s'K 



(4) 



4.2. Mass transfer rates 

There is no secure method of mass transfer rate determination for 
transient LMXBHs. One can obtain an estimate of this parameter 
by dividing the mass accreted during outburst by the recurrence 



The estimates of mass-transfer rates for leaky discs diff'er 
slightly from the ones given in Table 3 of Paper I, since in the 
present study we used in the equation for M^ax the lower limit 
of q as given in Table [T] instead of assuming a similar ^ = 0.1 
for all systems. For V406 Vul we estimated q by using the value 
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Figure 6. Masses of donor-stars in modelled population. Vertical 
lines show the ranges of M2 corresponding to the uncertainty in 
the determination of spectral types (Table [U. Heavy solid lines 
show "limiting" tracks as in Fig. |4] Thin solid lines are evolu- 
tionary tracks for 1 .0 and 1 . 1 M© donors with 4 accretors 
like in Fig. g] For XTE J1650-500 (at Po,-b=7.88 hr) there is 
only one determination of spectrum - K4V, but the same au- 
thors jOrosz et alj|2004 ') mention that the next best fits are G5V 
and K2III; for this reason we show the lower limit for M2 in this 
system. 



of the mass function and the mass of the secondary as obtained 
from its spectral type. 

The estimates of M,™x forXTE Jl 118+480, GRO J0422-H32, 
and GRS 1009-45 strongly suggest that the AML in short-period 
LMBHB might be defined by GWR only. 

H1705-25 (system 9) may, as noted in the previous subsec- 
tion, belong to the population in which donors fill Roche lobes 
when their central hydrogen abundance Xc is reduced by >50%. 

4.3. Masses of secondaries 

Figure |6] compares the ranges of the masses of donors in ob- 
served SXTs corresponding to the ranges of the estimates of 
their spectral types with the masses of donors in the model. The 
spec trum-mass scale is adopted after Schmidt-Kaler, as given in 
ICoxi(l2000). 

The situation with M2 is similar to that with Teff : our model 
population well covers the masses of the four shortest period sys- 
tems, but to explain longer period ones we need to apply tracks 
for more evolved systems that we did not include in our grid 
of tracks. Adding some AML due to MB also would improve 
the agreement with observations. Figure |6l like Fig.|4] suggests 
that the origin of LMBHBs with orbital periods of 10 - 12 hr 
may be associated with systems in which RLOF occurred when 
Xc < 0.35. 

5. Discussion and conclusion 

We have shown above that assuming the product of the stellar- 
envelope binding-energy parameter and the common-envelope 



expulsion efficiency parameter cce x /l=2, it is possible to re- 
produce, (within the uncertainty of observations) the number of 
LMBHBs in the Galaxy, the effective temperatures and masses 
of the donors in these systems (as inferred from the spectra of 
the latter) and their mass-transfer rates. This result is maintained 
for acE X /l=0.5 but further reduction of acE x /I to 0. 1 results in 
models whose parameters are not compatible with the currently 
available data on observed systems. Also we reiterate that (as 
found in Paper I) a substantial reduction of th e strength of mag- 
netic braking as compared to the "standard" dVerbunt & ZwaarJ 
1981) makes all calculated systems transient, in agreement with 
observations. 

The common envelope phase remains the most enigmatic 
phase of binary star evolution. As long as the processes of inter- 
action of the companion star with the envelope it is penetrating is 
not understood, the use of simple conservation-law based equa- 
tions will remain the necessary, albeit approximate, approach. 
However, the evolution of massive stars strongly depends also 
on stellar winds. The mass of the black-hole progenitor and the 
mass of its envelope, and the radius of the star that define the out- 
come of the common envelope stage, are interrelated via mass 
loss in the pre-common-envelope stage, which is not well con- 
strained. The survival of a binary in a supernova explosion de- 
pends on the mass-loss in the Wolf-Rayet star phase, the possi- 
ble kick imparted to the nascent black hole and the fraction of 
the mass of the exploding star that forms the black hole. None of 
these parameters are well constrained. 

The situation concerning mass-loss by massive stars is con- 
troversial. It became a recognized fact that stellar winds of both 
O-stars and WR- stars a re clumped (e.g., Owocki et al." 1988|; 
Smith 2007; Moffat et aL 1988; Hamann & Koesterke 1998) and 
empirical estimates of mass-loss rates that depend quadratically 
on the density have to be revised downward by a factor of 
several. In particular, downward revision of empirical M val- 
ues would bring them into agreeme nt with modern theoretical 
dVink et alj|200ll) rates for OB-stars dMokiem et alj|2007l; Ivin^ 



120071: IVink etal 



2007h. It is a lso claimed that the widely ac- 



cepted iNugis & LamersI (l2000l) rates for WR-stars, which al- 



ready are dumping-correc ted, have to be revised further down- 
ward dHamann et alj|2006l) . A decrease of M would mean more 
massive stellar envelopes and, generally, more energy would be 
needed for the ejection of common envelopes. This may be in- 
terpreted as a need for higher q-ce- 

On the other hand, based on results of extensive analysis 
of the ratios of blue to red supergiants, of Wolf-Rayet stars to 
O supergiants, of red supergiants to Wolf-Rayet stars and of 
the relative number o f Wolf-Rayet subtypes, WC to WN stars, 
Eldridge et al.l (1200 8') suggested that the total amount of mass 
lost by stars has to be increased. In the latter study both single 
and binary evolutionary mod els were considered and the mass- 
loss rates of lVink et all ( 1200 Ih and lNugis & LamersI d2000() were 
used. 

In such a controversial situation we may only claim, based 
on our results, that for the parameters of stellar evolution im- 
plemented in our evolutionary and population synthesis codes, 
which are consistent with state-of-the art stellar evolution the- 
ory, an agreement between the properties of the observed popu- 
lation of short-period LMBHBs and the model suggests a high 
efficiency of expulsion of common envelopes. 

Above, we have shown that the number and properties of ob- 
served short-period LMBHB may be explained by the model that 
assumes a strongly reduced strength of magnetic braking. But 
we rec all the existence of the alternative model of iMenou et al.l 
d 19991) which suggests that most LMBHB might have truncated 
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discs that are secularly in a cold and stable equilibrium with tran- 
siency due to random variabiUty in the properties of discs and/or 
mass transfer rates. Within this model most of the systems we 
modeled and classified as transients may appear as faint and sta- 
ble. At the moment there are no theoretical arguments against the 
Menou et al. model and observations that may serve as selection 
criterion between the two models still have not been defined. 
Thus, both models remain possible. 
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